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Abstract — Analytical solution performs a vital role in a wide variety of deep space exploration missions, especially the 
periodic solutions which were believed to be the unique avenue to solve three-body problem by Poincare. A,? the absence of a 
general solution for the problem, an approximate analytical solution of the circular restricted three-body problem is 
addressed by employing multiple scales method in conjunction with some analytical techniques. It is worthwhile to note that 
the presented solution in three-dimensional space contains two time scales: t and £t ( £ is a small, dimensionless 
parameter), which is significative to improve and perfect the known literature. 
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I. INTRODUCTION 

An existence of conditional periodic solutions of circular restricted three -body problem (CR3BP) was presented by Gao [1], 
However, because of lacking exact analytical solutions for the problem, a considerable amount of researchers were attracted 
to develop the approximate analytical solutions of the problem. In 1973, Farquhar and Kamel [2] developed an 
approximation method for this type of orbit. Richardson [3, 4] presented a third-order analytical solution about the collinear 
libration points of the CR3BP based on the method of successive approximations and a technique similar to the Lindstedt- 
Poincare method. Lu and Zhao [5] put forward a kind of improved third-order approximate analytical solution in 2009, which 
is more accurate than the classical analytic solution of Richardson. In addition, Nayfeh [6, 7] studied two types of resonance 
near the planar triangular libration points under a case of small amplitude. The planar 3:1 resonance with time scales t and 
2 

8 t was discussed when the planar restricted three-body problem was expanded to second-order terms in £ . Moreover, 
approximate analytical solutions of the planar 2:1 resonance with time scales t and £t was also investigated when the 
problem was expanded to third-order terms in £ . These creative results provide us with new motivation in studying the 
solutions of the three -body problem. 

In this paper, multiple scales method (see Nayfeh [8]) will be employed to construct a three-dimensional (3D) approximate 
analytical solution of a spatial CR3BP. The addressed 3D solution will be written in the following form 

U = £'U| (t, St) + £ 2 U 2 (t, St), (1) 

where U = (x, V, Z.) , U ( =(x i ,y i ,Z j ),i = 1, 2 , and £ is a small, dimensionless parameter. 

It is worthwhile to note that when the problem is expanded to third-order, the above presented solution in 3D space contains 
two time scales, which is significative to improve and perfect the known literature, in which the constructed solutions mainly 
possess three cases: a) The three -body problem discussed is a planar case, so the proposed approximate solution is a two- 
dimensional curve, b) The 3D approximate solution contains only one time scale when the spatial CR3BP was expanded to 
the third-order, c) The 3D approximate solution contains two time scales when the spatial CR3BP was expanded to the 
second-order. 


II. Construction of the 3D Approximate analytical Solution 

Note that the third-order approximate system with small amplitude of the CR3BP in rotating frame can be described as 
follows (the general third-order approximate system can be found at Koon et al. [9]) 


x-2y -(1 + 2c 2 )x = ~ £c 3 (2a 2 


2 2 

■y -z 


) + 2£c 4 (2x 2 -3y 2 -3z 2 ). 


(2a) 
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y + 2x + (c 2 - l)y = -3sc 3 xy - — sc 4 y (4x 2 - y 1 - z 2 ^, 
z + c 2 z = -3£c 3 xz--£c 4 z(4x 2 -y 2 -z 2 '), 

where 


(2b) 

(2c) 


1 . 



(±1 Tju+i-iy 


(i -ArT 
(i +r L ) n+l 


the upper sign is for L x and the lower one for L,, 


c„ = 


rl 


1 -JU + 


(1 +r L r 


for libration point L, , 


y L = n x , /?, is the angular velocity of relative movement between two primaries. 
The solution of equations (2a) ~ (2c) are assumed to possess the following form 


x = £x l (T () ,T { ) + £ 2 x 2 {T Q ,T l ), (3a) 

y = £y l (T 0 ,T l ) + £ 2 y 2 (T 0 ,T l \ (3b) 

z = £Z l (T Q ,T l ) + £ 2 z 2 (T Q ,T l ), (3c) 

where T (j = t, 7’ = £t . 

d d 2 2 5 

Then, time derivatives become — = D 0 + £D X H , — — = Dq + A£D 0 D X H , = , k = 0, 1, • • • . 

dr dr d T k 


Substituting equations (3) into equations (2a) ~ (2c) and equating the coefficients of£, £ , and £ to zero, leads to 


D 0 2 X[ - 2 D 0 y x - (1 + 2 c 2 )x x = 0, (4a) 

D 2 y x + 2 D q x x + (c 2 - l)Ti = 0, (4b) 

D hi +c 2 ^i =0- (4c) 

D 2 x 2 + 2D (l D l x l -2D Q y 2 —2D x y x -(1 + 2 c 2 )x 2 =0, (5a) 

D 2 y 2 + 2D l] D ] _y, + 2D 0 x 2 + 2D x x x + (c 2 -l)y 2 = 0, (5b) 

D o *2 + 2D <A r +c 2 z 2 = 0. (5c) 

2D q D\X 2 - 2 Dj 2 = fc 3 (2xf - yf -z 2 ), (6a) 
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2D Q D x y 2 + 2D x x 2 = -3c 3 x, , (6b) 

2D 0 D x z 2 = -3c 3 x x z v (6c) 

Let secular terms be equal to zero, the general solution of equations (4a) ~ (4c) and (5a) ~ (5c) read 


X 1 = p x COS(/10, (7a) 

y, = -Sp x sin(/t/), (7b) 

z, = p 2 cos(At)- (7c) 

x 2 = a 1 (T I )cos(yit + /?,(!))), (8a) 

y 2 =-Sa 1 (T l )sin(At + /? 1 (T 1 )), (8b) 

z 2 = a 2 (T x ) cos(/lf + /3 2 (T x )), (8c) 


where p x , p 2 and 8 are constants, A is the mould of the pure imaginary roots of equations (4), CC i and p j (1 = 1,2) 
satisfy 


a[{T x ) = 


8 (S-A) 


(( 2 + S 2 ) p\ - p\ ) cos(2 At) + ( 2 -5 2 }p{-p\ sin ( At + /?, (T x )) 


«,/?;oo= 


4(1 -AS) 
3 

8(8 -A) C 
3 


4(1 -AS) 


c 3 Spf sin(2/l0cos(/lf + /?,(7))), (9a) 

((2 + S 2 ) p\ - p ; ) cos(2A0 + ( 2 - S 2 ) p\ - p\ cos (At + /?, (T x )) 
c 3 8p\ sin(2 At) sin(At + /?,(7j)), (9b) 


3 

-2A [a 2 sin(/ir o + /? 2 ) + a 2 /? 2 cos(AT 0 + /? 2 )] = - — c 3 8p 1 p 2 [cos(2AT 0 ) + 1], (9c) 

Therefore, a third-order approximate analytical solution with two time scales to equations (2a) ~ (2c) can be represented as 

x = sp x cos(/li) + £ 2 a x (£t)cos(At + J3 x (£t)), (10a) 

y = s8p x sin(/lf) - £ 2 8a x (£t) sin(/li + /?, (£t)), (10b) 

z = £p 2 cos( At) + £ 2 a 2 (£t)cos(At + J3 2 (£t)), (10c) 


where Oi x and p x satisfy equations (9a) and (9b), Oi 2 and J8 2 satisfy equation (9c). 
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III. Conclusion 

Since the governing equations of the third body is time -dependent in inertial frame, which appear as a high-dimensional 
nonlinear autonomous system. According to the abundant known literature, this system seems unlikely to be solved via 
analytical approaches. For example, if we employ multiple scales method, then it is difficult solving the solutions to the 

2 

equations derived from equating coefficients of £~ to zero, not to mention the non-autonomous obtained from third -order 
terms in £ . However, these will be possible in rotating frame, where the equations are characterized by nonlinear 
autonomous system. In this frame, the method of multiple scales in conjunction with some analytical techniques is adopted to 
construct a 3D approximate analytical periodic solution for CR3BP. The solution was demonstrated with two different time 

3 

scales t and £t when the system was expanded to £ -order. This result improves and perfects the known literature. 
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